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We consider a short Josephson junction with a phase discontinuity k created, e.g., by a pair of tiny 
current injectors, at some point xo along the length of the junction. We derive the effective current- 
phase relation (CPR) for the system as a whole, i.e., reduce it to an effective point-like junction. 

From the effective CPR we obtain the ground state of the system and predict the dependence of its 
critical current on k. We show that in a large range of k values the effective junction behaves as a tpo 
Josephson junction, i.e., has a unique ground state phase po within each 27r interval. For k ~ 7r and 
xo near the middle of the junction one obtains a po ± p junction, i.e., the Josephson junction with 
degenerate ground state phase po ± p within each 2-7T interval. Further, in view of possible escape 
experiments especially in the quantum domain, we investigate the scaling of the energy barrier and 
eigenfrequency close to the critical currents and predict the behavior of the escape histogram width 
£t(k) in the regime of the macroscopic quantum tunneling. 

PACS numbers: 74.50.+r, 85.25.Cp 
Keywords: current-phase relation 


I. INTRODUCTION 

Recently a lot of attention is attracted to Josephson 
junctions (JJs) with an unconventional current-phase re¬ 
lation (CPR) 1,2 . In particular, <p 0 JJs 3-8 and ip JJs 9-15 
and their combinations 7 , proposed and/or demonstrated 
recently, show non-trivial physics 16 and have potential 
for applications in the classical 1 ' -19 and the quantum 
domains 20 similar to tv JJs. Here, tpo JJs are defined as 
JJs having a unique ground state phase (single Joseph¬ 
son energy minimum situated at) po ^ 0 within each 2tv 
phase interval, while p JJs (sometimes denoted also ±ip 
JJ) have a doubly degenerate ground state phase (double¬ 
well Josephson energy with minima at) ±ip within each 
2-7T interval. 

Currently, the classical properties of p JJs made of 
a short 0-7T JJ are understood rather well 12 ’ 13,19 . For 
example, p JJs have two critical currents J c _ and Ic+ 
corresponding to the escape of the phase from — p and 
+p wells. In our group we are starting investigation of 
quantum properties of such JJs. The first step in this di¬ 
rection could be an observation of the macroscopic quan¬ 
tum tunneling (MQT) of the phase 21-23 out of both —p 
and +p wells of the Josephson energy profile. For this 
purpose, one, usually, measures the phase escape statis¬ 
tics, by sweeping the bias current at a constant rate and 
measuring the exact value of the switching current many 
times. Assuming that in <p JJ at low temperatures (low 
damping) the initial state (—p or +p) is random 16 , the 
switching current histogram should have two peaks, each 
of them just below corresponding critical current I c ±. 
The widths cr(T) of each histogram peak usually (when 
the damping is small) decreases with decreasing temper¬ 
ature T. However, cr(T) is expected to saturate at some 


value (7 m i n for temperature below some T*. Such behav¬ 
ior is usually interpreted as a transition from the regime 
of the thermal activation of the phase over the barrier to 
the regime of the MQT of the phase through the barrier. 
However, it is necessary to show that the observed cr m i n is 
not related to the (noise in the) experimental setup and 
other trivial reasons. Usually, in such experiments one in¬ 
troduces some extra tuning parameter, e.g., a magnetic 
field, which allows to demonstrate that the setup is able 
to measure the histograms that are more narrow than 
cj m in- Simultaneously, for the MQT experiment with p 
JJ, it would be advantageous to have a tuning param¬ 
eter, which provides a continuous transition between tp 
(or ipo) JJ and a conventional 0 JJ, whose physics is well 
studied. 

For this purpose, we propose to use a short ID con¬ 
ventional 0 JJ, equipped with a pair of tiny current in¬ 
jectors. By sending a current I - m j from one injector to 
the other, we can create a k discontinuity of the Joseph¬ 
son phase at some point xq along the JJ, where injectors 
are attached 24-2 '. If tt = tv, the system is similar to 
a superconductor-insulator-ferromagnet-superconductor 
(SIFS) 0-7T JJ, which becomes a p JJ, if parameters 
are chosen correctly 28,29 . However, since n oc Ji n j is ad¬ 
justable, one can in situ tune the junction from a 0 JJ to 
a tp JJ and also study all the states in between. Tuning 
k one can also affect the widths cr m j n of the hystograms. 

The aim of this work is to develop a theoretical model 
for a short JJ with k discontinuity of the phase and to 
predict or interpret the results of MQT experiment such 
as the one outlined above. Namely, we derive an effective 
(averaged) CPR for a short JJ with a phase discontinuity 
k and obtain experimentally relevant quantities, such as 
the critical current or the escape histogram width as a 


2 


function of re. 

The paper is organized as follows. In Sec. II we intro¬ 
duce the model and present the averaged CPR and the 
averaged Josephson energy derived in details in appendix 
A. In Sec. Ill we obtain several experimentally relevant 
dependences such as the ground state phase, critical cur¬ 
rent and escape-related characteristics as functions of re. 
Sec. IV concludes the work. 


where 

Ujo(0) = - cos cos (6) - X 0 sin Q) sin(0); (7) 

U J2 (d) = sin2 (f ) cos ( 261 )- ( 8 ) 

III. RESULTS 


II. MODEL 


A. Comparison with the previous results 


We consider a short JJ of length 2w (w < 1 in units of 
Josephson length A j). The Josephson phase 4>(x) changes 
along the x coordinate (—w < x < +w). At x = Xq 
(— w < Xq < +w) there is a re discontinuity of the Joseph¬ 
son phase, created, e.g., by a pair of tiny (in theory in¬ 
finitesimal) current injectors 25 . The junction is biased 
by a uniform current density 7 (given in the units of the 
critical current density). Our aim is to derive an effective 
(averaged over the JJ length) current-phase relation for 
this system, i.e., "/(ip), where 

i\) = (</>( x)) = J <j>{x)dx, ( 1 ) 

is the average phase across the JJ. It is that is actually 
measured, if one considers the system described above as 
a black box with two electrodes. 

In appendix A we derive the averaged CPR of the sys¬ 
tem under question by using the perturbation theory up 
to the second order in w, treating w (the half-length of 
the JJ) as a small parameter. It is convenient to write 
the resulting averaged CPR as function of the phase 9 , 
which is related to the average phase ip across the JJ as 

0 = V<+|v o , (2) 


For re = 7 T we expect the result given by Eqs. (3), 
Eq. (4) and Eq. (5) to be similar to those previously 
obtained for a asymmetric 0-7T JJ 28 . To compare both 
results, first, we have to convert the phases. In Ref. 28 
the phase is continuous, while in our case it is has a 
discontinuity. Instead of the discontinuous phase we 

can introduce the continuous phase fi(x), which behaves 
exactly like <p(x), but without a re-jump, i.e., 

fi(x) for — w < x < Xq ^ 

n{x) + re for Xq < x < +w 

Then 

^<0 )=<M>+(l-*o)f. (10) 

It is (//) that is used in Ref. 28 (it is denoted as ip there). 
Rewriting our effective CPR in terms of (/z) and taking 
re = 7r, we obtain 

7 = V 0 sin(/z) - Qsin(2(/z)). (11) 

The quantities such as (j c ), L 0 and L, r from Ref. 28 can 
be expressed in terms of quantities used here as 


where Vo = Xq/w. The averaged CPR can be written as 
7(0) = 7o(0) + w 2 j 2 (0) + 0(w 4 ), (3) 

where 

7o(0) = cos (|) sin(0) - X 0 sin ^ cos(0), (4) 

is the 0 -th order result of the perturbation theory, the 
first order gives no correction, and 

72(0) = ^ sin2 (|) sin ( 2 0)> (5) 

is the second order correction in the framework of the 
perturbation theory, and Q = (w 2 / 6 )(l — Xq) 2 is intro¬ 
duced to make some formulas below more compact. The 
third order correction is zero and the terms 0(w 4 ) and 
smaller will be neglected. 

The effective Josephson energy of the system is an in¬ 
tegral of the effective CPR (3) and is given by 

Uj(9) = Ujo(6) + w 2 U J2 {9) + 0(w 4 ), ( 6 ) 


L 0 = w + x 0 = w{\ + X 0 ); 
L n = w - x 0 = w{ 1 - A' 0 ); 
Aq L-j, 


(jc) = 


Lq + L r 


= Vo. 


( 12 ) 

(13) 

(14) 


By substituting this into expression (18) of Ref. 28 and 
taking into account the definition of To, see Eq. (17) of 
Ref. 28, we arrive at the CPR (11) derived here. Thus, 
for re = 7 r the result of Ref. 28 is reproduced exactly. 
Note however that in Ref. 28 the small parameter is the 
deviation of the phase from its average value, while in 
our case the small parameter is w. Although, they are 
related (one expects small deviations for small w), this 
relation is not straightforward. 

In terms of variables used here 



(Vq 2 ~ l ) 2 


V 0 


(15) 


If r 0 < —1 then we have a ±ip JJ. This means that only 
for |A 0 | < u> 2 /3 one obtains a ip JJ at re = n. 
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FIG. 1. (Color online) The 7 c o(re) and (po(re) curves calculated using Eqs. (17) and (18) for (a) A'o = 0.5, (b) A'o = 0.2 and (c) 
X 0 = 0.05. 


B. Ground state phase and the critical current 

In the 0-th approximation the averaged CPR (4) can 
be rewritten as 

7 o = 7 co(«) sin (0 - 0 O ) = 7 co(«) sin(^ - ip 0 ), (16) 


vicinity of re = tt, see Fig. 2(b,c). To find the approximate 
analytical expression describing them we take re = tt + 5k 
(|5re| -C 1). Then we expand Eq. (19), up to the first 
order in 5k and solve it for 5k. Finally we obtain an 
approximate value of k for any given ground state phase 
0 gs , 


where 

7co(«) = \JXq sin 2 +cos 2 Q) (17) 

is the maximum supercurrent. The critical current mea¬ 
sured in experiment is ± 7 c o(re). 7 co(^) has maxima equal 
to 1 at k = 2irn ( n is any integer) and minima equal to 
Xq at k mod 27 t = 7 r, see Fig. 1. 

In Eq. (16) the ground state phase ipo = 0o — Xqk/2, 
where 


00 = arg 



+ i ■ Xq sin 



(18) 


and the function arg(z) returns the argument (phase an¬ 
gle) of a complex number 2 . Obviously, the CPR given 
by Eq. (16) corresponds to a ipo JJ 3 ~ 8 . Some examples of 
7 c o(k) and ipo(k) dependences are shown in Fig. 1. For 
large asymmetry Ao the modulation of 7 c o(re) is not as 
deep as for A'o —► 0. The phase shift ^q(k) changes from 
0 to 7 r(sgn(A 0 ) — X 0 ) as k changes from 0 to 27r. It is 
positive for Ao > 0 and negative for Ao < 0 . 

When A'o —> 0, the critical current given by the 0-th 
order formula (17) vanishes close to k = n and one has 
to take into account the next (second) order corrections 
given by Eq. (5). This happens for X 0 < w 2 / 3, see the 
discussion after Eq. (15). 

Next, we consider the second order approximation. 
The ground state phase 0 gs (re) is a solution of 7(0 gs ) = 0 
for this re, i.e., 

0 = cos sin(0 gs ) - A 0 sin Q j cos(0 gs ) + 

+ Qsin 2 sin(20 gs ). (19) 

This equation can be solved only numerically, see Fig. 2. 
It can be seen that multiple solutions 0 gs (re) appear in the 


re( 0 gs ) w 7r + 2 [2 Q cos( 0 gs ) - A'o cot( 0 gs )], (20) 

i.e., the inverse of the ground state phase 0 gs (re). This 
approximation is also shown in Fig. 2. One can see that 
approximation given by Eq. (20) is very good in the whole 
range of 0 < re < 2 - 7 t. Note, that the appearance of three 
solutions (two stable and one unstable) out of one near 
re = 7 T is a result of the competition of the cos-term with 
the cot-term in Eq. (20). From Eq. (20) one can fig¬ 
ure out that the multiple solutions appear for Xq < 2Q, 
i.e., |Ao| < w 2 / 3, which is in agreement with the dis¬ 
cussion after Eq. (15). We note that if Ao is so small 
it can be neglected in the definition of Q, so that for 
Ao ~ w 2 when the second order approximation becomes 
important, Q = w 2 / 6 . If |Ao| w 2 , one can omit the 
2Qcos(0) term ~ w 2 in Eq. (20) and end up practically 
with expression (18) for 0 o from the 0 -th order approxi¬ 
mation. The difference is that (18) gives only the stable 
branch, cf. the ground state phase shown in Fig. 1 (only 
stable branch) and Fig. 2 (both branches). 

From Eq. (20) one can find the range of re where the 
double ground state exists, i.e., the points Q\ and 02 in 
Fig. 2 where dre(0 gs )/d0 gs = 0. We obtain that 



which lay symmetrically with respect to 0 = 7 r/ 2 . It fol¬ 
lows from Eqs. (21) that the bifurcation point, where 
the 0 (re) curve switches from one stable to three so¬ 
lutions (two stable and one unstable), corresponds to 
0i = 02 = 7 t/2 , i.e., at A 0 = 2Q « w 2 / 3, which is 
again in agreement with the result obtained directly from 
Eq. (20). The range of re around re = 7 r where two sta¬ 
ble solutions exist, is found by substituting Eq. (21) into 
Eq. (20). 























4 



FIG. 2. (Color online) The ground state phase # gs (re). Comparison of the approximate dependence given by Eq. (20) (black 
line) with the exact dependence calculated by numerically solving 7 ( 0 ) = 0, see Eq. (3), for each value of re (symbols), (a) 
Xo — 0.10, (b) Xo = 0.04, (c) A'o = 0.04 and the region close to re = n zoomed. In all plots the regions with a positive slope 
(green) correspond to a stable solution (energy minimum, where U"(9) = 7 '(d) > 0), while the regions with a negative slope 
(red) correspond to an unstable one (energy maximum, where U"(8) = q'( 6 ) < 0). 


Looking at the Fig. 2 one sees that the ground state 
phase 9 gs has the values symmetrically placed around 
9 = (7t/ 2) sgn(Xo) only at re = n. From Eq. (20) they 
are given by 


6» gs i(7r) = arcsin 


*0 


2 Q 


and 0 gS 2 ( 7 r) = it — arcsin 


Xn 


2QJ 

( 22 ) 

For re 77 the symmetry is brocken because the corre¬ 
sponding double-well potential Eq. (6) becomes asym¬ 
metric (one well is deeper than the other) relative to 
9 = ( 77 / 2 ) sgn(Xo). The real (measurable) ground state 
phases are given by ip gs = 9 gs — Xqk/2 and at re = 77 are 
symmetric with respect to the phase (7r/2)[sgn(Xo) — Xo]. 
Recalling that for doubly generate state to occur one 
need very small |Xo| < w 2 , the shift from ( 77 / 2 ) sgn(Xo) 
is small. Thus, such a JJ can be called a <po ± ip 
JJ, where, at re = tt ipo = (77/2)[sgn(X 0 ) — X 0 ] and 
ip = 77/2 — arcsin(Xo/2Q). 

To find the critical current (s) in the second order ap¬ 
proximation for each re we search for an extremum of 
7 (0, re) with respect to 9. It takes place at 9 = 6 C for 
which 


7 '(0c, k) = X 0 sin Q j sin(0 c ) + cos cos(0 c ) + 

+ 2 Q sin 2 ^ [2 cos 2 ( 6 » c ) - 1] = 0. (23) 

Here and below the prime denotes d/d9. This equation 
can be solved for 9 C only numerically to find several (up 
to 4) 9 C for each value of re. Then we substitute each of 
these 9 C into Eq. (3) to find y c (re) = 7 (9 C , re). The result 
is presented in Fig. 3. The global behavior is defined 
mainly by 70 , i.e., j c o. However, near re = 77, where 7 c o 
vanishes, 72 results in a bistability and in the formation 
of a <0-like intersection of the branches. Such behavior 
is typical for a p JJ made of 0 and 77 parts 7,28,29 . 

Similar to the case of the ground state phase, one can 
find an approximate expression for 7 c (re) near re = 77. By 


substituting re = 77 + Sk (|dre| -C 1) into Eq. (23), Taylor- 
expanding it up to terms 0(5n ), and expressing 8k, we 
obtain the critical value of re corresponding to y c for given 

9. 


K c(0) =TT + 8KKiTT + 2 


Xo sin( 0 ) + 2 Q cos(2 9) 
cos (9) 


(24) 


To calculate the critical current, we substitute Eq. (24) 
into Eq. (3), which was preliminary expanded near re = 77 
up to O(dre). We obtain 


7 c P ( 0 ) ~ -^ocos( 6 >) + Qsin( 20 ) — 

X 0 sin 2 (0) + 2 Q sin (9) cos(2 9) 
cos (9) 


(25) 


By sweeping 9 in the range —77 ... 77, we can now calculate 
k c (0) and 7 c P ( 0 ) and make a parametric plot of 7 c P ( 0 ) vs. 
re c (0), see Fig. 3. The agreement with direct numerical 
calculations near re = 77 is excellent, see Fig. 3(b,c). The 
deviations become noticeable as re mod 277 approaches 0 
or 277, see Fig. 3(a). 

This approximate analytical expression for 7 C allows 
us to calculate some key features in the 7c(re) plot. For 
example, one can find out the value of d<> (and re<>), for 
which the branches meet each other, see points A and 
B in Fig. 3. The analysis of the re c (0) dependence (24) 
shows that this happens when dn c /d9 = 0. Differentiat¬ 
ing Eq. (24) we obtain the following equation for 9 <>. 


4<5sin 3 (0^) — 6 Qsin( 6 | o)+X 0 = 0. (26) 


This cubic equation with respect to sin(d<>) has only one 
suitable root, which (after some lengthy algebra) can be 
expressed as 


sin(0o)=y= c °s(f)-vising) , (27) 


where \ ma Y be explicitely written as 


X = arg 


—Xq + i ■ 



(28) 
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FIG. 3. (Color online) Examples of the 7 c (k) dependence for w = 0.5. (a) A'o = 0.04, global behavior of 7 c (k); (b) Xq = 0.04, 
zoom of the region of interest near k, = n; (c) A'o = 0 . 02 , zoom of the region near k — n. Thick (red) lines/symbols show 7 C 
obtained by directly solving Eq. (23) numerically to find all 9 C and then calculating -y c from Eq. (3). Thinner black lines/symbols 
correspond to the approximation given by Eq. (25). Gray dashed lines show 7 c q(k), see Eq. (17), for the same parameters. 


It turns out that for typical parameters corresponding to 
a bistable case (small |Ao| < w 2 /3 ~ 2 Q), the value of x 
changes from 7r/4 to 37r/4. Then by expanding Eq. (28) 
as 


In essence this is a CPR. The critical current is reached 
for 9 = 9 C , when d r y/d9 = 0, i.e., 

U'j(e c ) = 0. (33) 


7 r V2Xo 7 r 3 Xo 

2 + ^r~2 + 7i^’ 


(29) 


we get sin(0<>) « X 0 /(6 Q) = X 0 /w 2 < 1/3. In fact, 
this limit of small sin(6^) corresponds to neglecting the 
sin 3 (0<>) term in Eq. (26), so that one obtains sin(0<>) « 
Xq/(6Q) right away from Eq. (26). After finding sin(0<>), 
the value of k<> can be found as K c (0q) from Eq. (24). Us¬ 
ing the approximation (29), i.e., in the worst case neglect¬ 
ing 4sin 3 (0<>) = 4/27 in comparison with 6sin(0<>) = 2 
(accuracy ~ 8%) in Eq. (26), we can write 


K<> 


7r ± 


^o 2 


(30) 


C. Energy barrier 

We consider the thermal escape or the quantum tun¬ 
neling of the phase 9 out of the potential well, when the 
bias current 7 —> Since our model reduces the 

system to an effective point-like JJ, for calculation of the 
escape rate T one can use standard thermal or quasi- 
classical quantum formulas. In these formulas, the key 
parameters are the barrier height AU(y) and the eigen- 
frequency wo ( 7 ). The aim of this section is to obtain the 
expressions for them. 

In general, we proceed as follows. Given the Josephson 
energy profile Uj(9 ), the total potential energy of the 
biased JJ can be written as a tilted potential 


From here one can, in principle, find (one or more) values 
of 9 C . Imagine that we have found all values of 9 C . Then, 
the value of the critical current y c is found from Eq. (32), 
i.e., 


7 c = U'j{0 c ). (34) 

Now we assume that the value of 7 is slightly undercrit- 
ical, ie., 7 = 7 C (1 — J 7 ), where 0 < J 7 -C 1 for any sign 
of 7 c . We are interested to expand the potential U(9) in 
the vicinity of the bending point 9 C . We, therefore, write 
9 = 9 C + 59 (|<50| <C 1) and substitute this into Eq. (32) 
and Taylor-expand up to 0(69 2 ). We get 

U'AO c) + Uj(@c)59 + \u'j\9 c )59 2 = 7c (l - <J 7 ). (35) 

Here the first terms in the l.h.s. and the r.h.s. cancel 
because of Eq. (34), the second term in the l.h.s. vanishes 
because of Eq. (33). As a result we obtain new static 
solutions for an undercritical 7 shifted from 9 C by 


to the positive or negative direction. One of them is 
stable and corresponds to the minimum of U(9), another, 
unstable one, corresponds to the maximum of U(9). The 
energy barrier 

AU = \U{9 c + 69)-U{9 c -59)\. (37) 


U(9) = Uj(0)- 7 0 . 
The static solution(s) correspond(s) to 

U'AO) -7 = 0, 


(31) After expanding up to 0{59 3 ), we see that the terms 
0(59°) and 0(59 2 ) cancel, and we obtain 


(32) A U = \2UAlc)59 - 2^69 + -U'f(9 c )69 3 \. (38) 
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Using Eq. (34), the definition of 7 = 7 C (1 — £ 7 ) and the 
expression (36), we finally obtain 


4\/2 l7c^7l 3/2 

3 


4V2 

— 



(39) 


Now let us apply this general result to our system. 

If Xq w 2 /3 or if k is far away from k = 7 r, then 

we can use only the 0-th order term in Eq. (3) and 
in Eq. ( 6 ). In this limit the CPR is sinusoidal, see 
Eq. (16). Although it is shifted by 9 0 (<po), it is irrel¬ 
evant for calculation of the escape barrier and eigenfre- 
quency. Thus, the system behaves as a conventional JJ 
with sinusoidal CPR, critical current 7 c o(/i) and Joseph- 
son energy Uj = 7 c o(« 0 [l - cos(0 - 0 O )]. Thus, 0 C = 7 r/ 2 , 
U"'(6 C ) = — 7 co and we obtain the usual approxima¬ 
tion for energy barrier and eigenfrequency in the limit 
7 ->• ±7co: 


7 r the positive 7 c (k) branch crosses zero, see Fig. 3(b), 
so that we see that the prefactor also vanishes at this 
point, see Fig. 4(b). Then, at somewhat larger k both 
mentioned branches join, see Fig. 3(b). At this point the 
prefactor diverges, see Fig. 4(b). The other two branches 
in Fig. 4(b) show similar behavior. 


D. Eigenfrequency 

In general, the eigenfrequency of phase oscillations 
around one of the static solution 0 C ± SO, see Eq. (36), is 
given by 

uj 2 = U"{0 C ± 50) = U"(0 C ) ± U"'(0 c )S0. (44) 

The first term vanishes because of Eq. (33). So, using 
Eq. (36) we get 


AU (£ 7 , k) = -^7 cO («0 Ay 3 / 2 . (40) 

In the case when the second order correction is im¬ 
portant, i.e., Xq < ur/3 and k « ir, we use the same 
approach, but again, like in the section about critical 
current and ground state phase, we approximate for 
k = 7 r + 5k (|5k;| -C 1). In this case the energy is given 

by 


to 2 = V| 27 ct/"'( 0 C )| 5^ 2 . (45) 

In the 0-th approximation 0 C = n/2, U"'(0 C ) = — 7 c o 
and we arrive to the well-known result 

u, 2 = lc0 V26^ 2 . (46) 

In the second order approximation we again sweep 0 
and make an implicit plot of the Wq prefactor 


& K, G 

U (0, Sk) = — cos(0) — A 0 sin( 0 ) — cos( 20 ). (41) 

From here 

U'"(6, k) = — sin(0) + Xo cos(0) — 4 Q sin(20). (42) 

We have to take 0 = 0 C in Eq. (42) and substitute 
this into Eq. (39). The dependence of AU on Ay is ob¬ 
vious from Eq. (39), so our aim is to see how prefactor 
in Eq. (39) depends on k (or 5k). Since, 0 C for each k 
can be found only numerically, we, as in the previous sec¬ 
tions, sweep 0 from — 7 r to + 7 r and find the corresponding 
k c (0) from Eq. (24) and then calculate U"'(0 , k c (0)) from 
Eq. (42). Then we make a parametric plot of the energy 
barrier prefactor 


-Pac/(0) = 


A U{0,k c {0)) 

,5y3/2 


4V2 

— 


7 c ( 0 , K e ( 0)) 3 

C7'"(0,K c (0)) 


, (43) 


as a function of k c ( 0), see Fig. 4. The global behavior 
is given by the 0-th order approximation, see prefactor 
in Eq. (40). The second order approximation, where we 
expanded all expressions near k = ir, works well near k = 
7 r, but deviates substantially from the real solution given 
by the 0-th order approximation when (k mod 27r) —> 
0 or —> 27r. In Fig. 4(b), as k increases, one sees two 
branches, given by Eq. (43). One of them corresponds 
to the negative critical current branch, cf. Fig. 3(b), 
another to the positive one. At k slightly larger than 


PZ o (0) = ul/5^ 12 = V|2 7 c (0,« c (0))J7"'(0,« c (0))|, 

(47) 

as a function of « c (0) given by Eq. (24). The behavior of 
the eigenfrequency prefactor is shown in Fig. 5. Similar 
to the energy barrier prefactor, the eigenfrequency pref¬ 
actor given by Eq. (47) describes the multiple solutions 
near k = it well. However, the 0-th approximation is 
better outside this vicinity. 


E. Escape histogram width in the MQT regime 

The dependences AU(5"/,k) and ui o(Ay, k) allow to di¬ 
rectly calculate not only the escape rate, but also the 
width g of the escape histogram as a function of k. This 
g(k) dependence can be directly compared with the ex¬ 
perimentally measured one. For the sake of simplicity we 
limit ourselves to the case of MQT, so that the temper¬ 
ature is excluded. The approximate, but rather precise, 
formula for the histogram width (dispersion) cr was de¬ 
rived by Garg 30 in the general case of a particle in a 
tilted potential. For MQT regime the Garg 30 expression 
reduces to 


(48) 

where we have omitted ln-terms that are much weaker 
than power-terms. We use a oc sign as we are interested 
not in the width itself but in its scaling as a function 


( 7 OC 


Pu, 


Pau 
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FIG. 4. The energy barrier prefactor Pau given by Eq. (43) vs. k c (9) (black) and the one calculated using the 0-th order 
approximation Eq. (40) (gray), (a) shows global behavior in the interval 0 < k < 27t, while (b) shows the zoom of the area 
close to k = 7r, where multiple solutions appear. Parameters are: w = 0.5, Xo = 0.04. 



FIG. 5. The eigenfrequency prefactor P% 0 given by Eq. (47) vs. k c (6) given by Eq. (24) (black) and the one calculated using 
the 0-th order approximation Eq. (46) (gray), (a) shows global behavior in the interval 0 < k < 2-7T, while (b) shows the zoom 
of the area close to k = 7r, where multiple solutions appear. Parameters are: w = 0.5, Xq = 0.04. 



FIG. 6. The width <t(k) of the escape current histogram in MQT regime calculated using the 0-th order approximation Eq. (50) 
(gray line) and the second order approximation linearized near k = n, see Eq. (49). (a) shows the global behavior in the interval 
0 < k < 2n, while (b) shows the zoom of the area close to k = ir, where multiple solutions appear. Parameters are: w — 0.5, 
X 0 = 0.04. 



























of k. This a is a dispersion of S 7 defined above, i.e., 
it assumes that the critical current is equal to 1. If the 
critical current is equal to 7 C , the sigma (measured in the 
same units as 7 c ) is 


cr(re) 


Pui 0 (k) 

_Pau(k,)_ 


4 

5 


l7c(«0|. 


(49) 


When Xq w 2 / 3, or k mod 2n is not very close to 7 r, 
we can use the 0-th order approximation. In this case, 
by substituting the prefactors from Eq. (40) and Eq. (46) 
into Eq. (49), we get 


desired xq — a small technological shift can drastically 
change the junction. 

However even a nominally symmetric junction (Xq = 
0 ), due to a tiny technological misalignment can get 
Xq ^ 0. As a result an experimental |7 c (k)| curve ex¬ 
hibits asymmetric minima for positive and negative bias 
current, as in Fig. 3. Also, in experiment it should be 
easy to conclude whether the asymmetry |Xo| is smaller 
than w 2 /3 (and we deal with a ±</? JJ) or larger (and 
we deal with single state ipo JJ). In the former case the 
7 C (k) dependence should have a cusp-like minimum with 
branch crossing, while in the latter case the minimum 
will be smooth. 


(J 0 (k) oc 7^0 5 («0 ■ 


(50) 


IV. CONCLUSIONS 


This dependence is shown by the gray line in Fig. 6 . 

For small Xq < w 2 /3 and k mod 2n in the vicinity 
of 7 r, we have to use the second order formulas. Again 
by substituting prefactors from Eq. (43) and Eq. (47) 
into Eq. (49) we obtain a(9), which we plot vs. n c {6) 
as a parametric plot, see Fig. 6 . One can see that a 
vanishes at the bifurcation point where the two branches 
join. Note also that at the points k = k z where 7 c (k) 
vanishes (crosses zero), i.e. 7 C oc n z + <5k is linear, both 
Pau oc <5k 3 / 2 and P Uo oc (5k 1 / 4 vanish, however er oc 
const does not have zero or any other peculiarity at these 
points, as can be seen from Eq. (49). 


F. Experimental relevance 

On one hand the range of |Ao| < w 2 /3 (|xo| < w 3 / 3), 
required to create a (/? JJ near k = 7 r, is very tiny. This 
was already pointed out in the previous works' -11,13,28,29 . 
This makes it very difficult to controllably fabricate the 


We have derived an effective model, which describes 
a short JJ with a phase discontinuity k at an arbitrary 
point xq along its length. This model reduces the system 
considered here to a point-like JJ with an unconventional 
CPR. One can relatively easy obtain all desired charac¬ 
teristics of such a point-like JJ. For example, we analyzed 
the ground state and found that close to k = 7r one ob¬ 
tains a <fQ±ip JJ, while far from this point it is tpo JJ. We 
also calculated the dependence of the critical current of 
such a JJ as a function of k and found multiple branches 
close to k = 7T, corresponding to ipo ± ip states. Further, 
we have calculated the behavior of the energy barrier and 
eigenfrequency close to the critical current, which allow 
to make estimations of the width cr of the switching cur¬ 
rent histogram in the regime of macroscopic quantum 
tunneling. 

Appendix A: Derivation of the averaged CPR 

We use the following phase ansatz 





A l (x - Xq) + \B l (x - Xq) 2 + \C L (x - Xq) 3 + \D L (X - XqY , 
A R (x - * 0 ) + \B R (x - Xq) 2 + \C R (x - (To ) 3 + \P>r{x - Xq ) A , 


— W < X < Xq 
Xq < X < +W ’ 


(Al) 


which corresponds to the phase discontinuity at x = xq 
and a Taylor expansion of the phase 4>(x) in the left and 
right region (subscripts L and R) relative to the discon¬ 
tinuity Xq. 

In statics, the phase cj>(x) should satisfy the Ferrel- 
Prange equation 

7 = sin (4>) - </>", (A2) 

subject to boundary conditions at the edges x = ±w, cor¬ 
responding to zero applied magnetic field, and the held 


continuity at x = Xq: 


<p'(—w) = 0; (A3a) 

cj)'(+w) = 0; (A3b) 

<!)'(xq - 0) = 4>'(xq + 0). (A3c) 


The prime denotes d/dx. Below, we use the junction 
half-length w as a small parameter and develop a pertur¬ 
bation theory with respect to w. 







9 


1. O-th approximation in w 

We substitute the ansatz Eq. (Al) into the Ferrel- 
Prange Eq. (A2). After calculating <j )", we would like 
to expand sin(</>) with respect to the small parameter w. 
The key point is to make this expansion correctly. For 
this, we transform the argument of the sine function to 
explicitly pull out w from all terms. Namely, we define 
that Al,r and Bl,r from ansatz Eq. (Al) depend on w 
as 

x — xq = A Xw\ (A4) 

A l ,r = a L i,Riw, (A5) 

Bl,r = bLo,Ro- (A6) 

Here and below the subscripts separated by a comma 
mean that it is actually two equations: one is obtained 
by taking the first subscript in the whole equation, the 
second equation is obtained by taking the second sub¬ 
script in the whole equation. The higher order C and 
D-terms from Eq. (Al) are not relevant in O-th approx¬ 
imation. Initially, the scaling of A and B with w is ac¬ 
tually not obvious, but later we will see that the scaling 
given by Eqs. (A5) and (A6) is consistent. After the 
above substitution we expand sin(...) in Eq. (A2) rela¬ 
tive to w, keeping only constant terms (neglecting 0(w) 
and smaller). We arrive at the following expression(s). 

bw,Ro = sin (^o O' 0 ’ (A7) 

where 7 = 7o in our O-th approximation. From the 
Eqs. (A7) it is obvious that B scales ~ w° as written 
in Eq. (A6). From the boundary conditions Eq. (A3) we 
have 


Oil = (Alo + l)&L0i 

(A8a) 

a-Rl = (-^0 — 1)&lo; 

(A8b) 

Oli = Ofll) 

(A8c) 


where Xq = xq/w. It is Eqs. (A8a) and (A8b) where 
it becomes obvious that Alr ~ w, as it was correctly 
written in Eq. (A5), otherwise the l.lr.s. and the r.h.s. 

I 


bw (, « 

a L 3 = - 7 j- cos (^ 0 - - 

bRO (, , «A 

aR3 = - 7j- cos (^>0 + 2 ) 

®L3 = UR3- 


would have different orders in w. By substituting 6 lo,ao 
from Eqs. (A7) into Eqs. (A8a) and (A8b) and then a^i 
and clri from Eqs. (A8a) and (A8b) into Eq. (A8c), we 
finally get the current-phase relation 

7o(</>o) = cos sin(0 o ) - X 0 sin cos(0 o )- (A9) 

2. 2-nd order approximation 

In the next order (~ w 2 ) approximation we use all the 
terms in ansatz (Al) to substitute into the Eq. (A2). Af¬ 
ter calculating (f>" 1 we explicitly extract w from all terms 
using the following substitutions 


Al,R = aLl,RlW + aL3,R3W 3 ', 

(AlOa) 

Bl,R = bw,R0 + bL2,R2W 2 -, 

(AlOb) 

Cl,R = CLl,RlW, 

(AlOc) 

7 = 70 +w 2 72 . 

(AlOd) 


Here cili,ri, &lo,.ro an d 70 are from the O-th order ap¬ 
proximation and a,L3,R3 and 6l 2,A2 are the next order cor¬ 
rections. Other powers, e.g ., a,L2,R2 = = 0. After 

the above substitution we expand sin(...) in Eq. (A2) 
relative to w, keeping only terms ~ 0(w 2 ) and larger 
(neglecting 0(w 3 ) and smaller). We arrive at the second 
order polynomial in AX equal to zero. Obviously, it can 
be equal to zero for any AX only if each coefficient in 
front of AX 2 , AX and constant are all equal to zero. 
We, thus, obtain 

AX° : 6l2,H2 = —72; (Alla) 

AX 1 : c L i,ri = -a L 1>R1 cos (</> 0 T ; (Allb) 

AX 2 : D l ^ r = —b L Q R Q cos =F ■ (Allc) 

As a next step we substitute the ansatz (Al) into the 
boundary conditions (A3). Then we substitute the defi¬ 
nitions (A10), cancel the terms from O-th order approxi¬ 
mation (if any), substitute Cli.ri, D rr and bL2,R2 from 
Eqs. (All) and obtain 


(X 0 + 1) 3 -72(Xo + 1); 

(Al2a) 

(Xo-1) 3 -7 2 (Xo-1); 

(A12b) 


(Al2c) 


By substituting Eqs. (Al2a) and (A12b) into Eq. (A12c) we finally obtain 


1 



K\ 

72 = g cos 1 

(*+ 2 ) 

1 b R0 - - cos ^ 

l<M 

0 

e 
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Finally, we substitute o,ro from Eq. (A7) (0-th approx¬ 
imation) and obtain the final expression for the second 
order correction to the current 

72(</>o) = j;(l - ^o) sin2 (tj) sin(2<^ 0 ) + 

+ ^o(l - -^o) sin («) [l + cos(2^ 0 )], (A14) 


3. From <i>o to the average phase ip 

Up to now both the 0-th order CPR Eq. (A9) and the 
second order CPR Eq. (A9) are given as functions of the 
(po, while our aim is to express those as a function of 
the average phase ip. To find ip, we substitute ansatz 
Eq. (Al) into Eq. (1), integrate on each interval, sub¬ 
stitute the definitions of A , B , C, D from Eqs. (A10), 
then substitute expressions for a, 6 , c, D and 70 from 
Eqs. (A7), (All), (A12) and (A9). Then we keep only 
the terms ~ 0(w 2 ) and larger, after some simplifications 
arrive at 

= <po-^X 0 + w 2 ‘^X 0 (l-X$)sm Q) cos{(po). (A15) 

Our aim is to invert this expression, i.e., to express <pa{ip) 
to substitute to CPR and obtain 7 (ip). We again act fol¬ 
lowing the perturbation theory with respect to the small 
parameter w. In the 0-th approximation 

<^ 0) = ip + ^x 0 = e. (Aie) 

Here we introduced the angle 9 , which makes expressions 
more compact. 

In the next (second) approximation (po = cp g 0 ' +w 2 (P < 'q\ 
By substituting this into Eq. (A15) and expanding up to 

0(w 2 ), we obtain cp^ and therefore 

<Po(4>) = 9-w 2 ^X 0 { 1 - Xq) sin Q) cos (6). (A17) 

Finally, we substitute this into expressions (A9) and 
(A14), expand up to 0{w 2 ) and after some simplifica¬ 
tions obtain 

7 (9) = cos (0 sin(fl) - X 0 sin Q) cos( 6 ») + 

+ Q sin 2 sin(20), (A18) 

where, for the sake of brevity, we have introduced 

q= l w ^i-X 2 ) 2 . (A19) 
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